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ABSTRACT 

The stellar scintillation is one of the fundamental limitation to the precision of ground- 
based photometry. The paper examines the problem of correlation of the scintillation 
of two close stars at the focus of a large telescope. The derived correlation functions 
were applied to data of the long-term study of the optical turbulence (OT) in the 
Northern Caucasus with MASS (Multi- Aperture Scintillation Sensor) instrument to 
predict the angular correlation of the scintillation at the Sternberg institute 2.5 m 
telescope currently in construction. A median angular radius of the correlation as 
large as 20 arcsec was found for the case of Kolmogorov OT. On the basis of the 
obtained relations we also analyze the correlation impact in ensemble photometry and 
conjugate plane photometry. It is shown that a reduction of the scintillation noise up 
to 8 times can be achieved when using a crowded ensemble of comparison stars. The 
calculation of the angular correlation can be repeated for any large telescope at the 
site where the OT vertical profiles are known. 
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1 INTRODUCTION 

Many astronomical problems require high-precision photo- 
metric data. Two clear examples are the determination of 
parameters of ext raterrestrial planetary systems from eclips- 
ing light curves (lEverett fc Howellll200ll ; IWinn et al.ll2009l ; 
ISouthworth et all l2009l h~ and the study of low-amplitude 
stellar pulsations, reflecting inner structure of stars 
( Heaslev et al.fl996l ; lAerts. Christensen-Dalsgaard fc Kurtj 
l2010l ). The success of the CoRoT and Kepler space missions 
clearly demonstrates the need for further increase of the ac- 
curacy of astronomical photometry. 

Despite the superiority of space projects, observations 
with ground-based telescopes also have some advantages. 
They are more accessible, the number of the telescopes is 
large and their service life is tens times larger than for space 
missions. The larger apertures of these telescopes push the 
intrinsic limit of the accuracy, caused by photon noise, to- 
wards fainter objects or improved temporal resolution. 

The main limitation to obtain high-precision photomet- 
ric results from the ground is due to the Earth's atmo- 
sphere, in particular the atmospheric scintillation of stars. 
The stellar scintillation is a random temporal fluctuations 
of the radiation entering the telescope aperture and is re- 
lated to the amplitude distortions of the light wave having 
passed through the turbulent atmosphere. Photometric er- 
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rors induced by the scintillation (scintillation noise) have re- 
peate dly been studied (see, e.g.. lYoundll969l ; lDravins et al.l 
Il997ft , 

In comparison to the photon noise, the scintillation 
noise does not depend on the brightness of the stars and 
thus cannot be reduced by choosing brighter objects. For 
typical exposure times (few seconds and longer), the scin- 
tillation noise decreases more slowly with telescope diam- 
eter than the photon noise and thus becomes to dominate 
on large telescopes. Consequently exactly t he scintillation 
noise sets a limit to the photometric accuracy l|Heaslev et al.l 

1996) . This problem becomes even more relevant for expo- 
sures shorter than 1 s. 

In contrast to the flux variations due to changes in at- 
mospheric transparency, the stellar scintillatio n is character- 
ized by a lower spatial (angular) correlation (|Dravins et al.l 

1997) . That leads to an increase of the errors in differen- 
tial photometry when using a distant comparison star. It 
is assumed that the angular correlation for the scintillation 
is typically a few arcseconds but the special case of large 
telescopes remains virtually unexplored. 

This paper is a theoretical study of the angular correla- 
tion of the stellar scintillation on the basis of researches on 
the scintillation as a tool for probing OT above astronomical 
sites. We first consider the link between OT parameters and 
photometric scintillation noise. Then the formulas for the 
correlation at large telescopes in two limiting regimes, short 
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(much less than 1 s) and long (greater than 1 s) exposures, 
are derived and analyzed. 

In Section [3] the results are applied to the OT vertical 
distribution above M t Shatdzhatmaz in the North Caucasus 
l|Kornilov et al-lfeoiOT ) where the Sternberg Astronomical In- 
stitute (SAI) 2.5 m diameter telescope is to be installed. The 
angular correlation functions and typical angles of the cor- 
relation are calculated. Finally, the results are applied to 
two methods aimed to reduce the scintillation noise, the re- 
cently proposed pho tometry in a conjugated entrance pupil 
dOsborn et alj201Ch and the more traditional ensemble pho- 
tometrv i|Girliland et all 1 1993ft . 



2 THEORY 

2.1 Scintillation noise in photometry 

In differential photometry, the light flux from the science 
object 1 is compared to the flux from another reference ob- 
ject 2 in order to remove the uncertainties caused by in- 
strumental factors and by variations in atmospheric trans- 
parency. The comparison can be performed as the ratio of 
I1/I2 but more often the difference between stellar magni- 
tudes: Am = 1.086 \nh/h. 

Several noise sources of different nature affect these 
measurements. Depending on its origin, the noise enters in 
the measured signal additively (detector noise, sky back- 
ground noise), multiplicatively (atmospheric transparency 
fluctuations and scintillation noise) or in a more complex 
way (photon noise). Hereafter, we consider only one compo- 
nent, the scintillation noise SI which is proportional to the 
mean light intensity or its equivalent in magnitudes, the ad- 
ditive noise 5m. By definition, the mathematical expectation 
(SI) = and (8m) = 0. 

Generally accepted that the power of the stellar scin- 
tillation is measured with the scintillation index which is 
defined as s 2 = (SI 2 )/(I) 2 or s 2 = {(In I - ln(/)) 2 ). For 
small fluctuations of flux, these definitions are equivalent. 

Therefore, we can present the noise variance (scintilla- 
tion noise power) in the form as erf = (-O^ 2 f° r intensities, 
or cr 2 = 1.179 s 2 for magnitudes. The normalised variance 
of the ratio o = I1/I2 and the variance of the magnitudes 
difference o\ m are simply expressed in the terms of the scin- 
tillation indices: 



„2 



Si +S2 - 2ci, 2 , 



(1) 



(2) 



o 2 1.179 

where ci,2 is the cross-index, defined naturally as 
_ (8h 8h) 

Cl ' 2 ~ (h)(h) ' 

When two objects are measured simultaneously in the same 
field of view, we can assume that their scintillation is equal: 
s'{ — s 2 , and the cross- index depends only on the angular 
distance 8 between the objects. Thus, the previous formula 
(fTJ can be written as s d = 2 s — 2 c(9) . 

For uncorrelated signals, s 2 is twice as large as the scin- 
tillation index, i.e. the scintillation noise will be doubled 
after the comparison. On the contrary, in the case of a sig- 
nificant correlation, the scintillation noise is reduced. This 
is a very important fact which, therefore, has been exam- 
ined iri_botli_y2«ore_tical (|Rvan fc Sandlerlll99Sft and practi- 
cal (|Dravins et al.lll997ft aspects. 



It seems reasonable to define the angular correlation 
radius 8\ as the angle at which sj| becomes equal to s 2 , i.e. 
when there is no amplification of the scintillation noise after 
to comparison. Obviously then c(9\) = s 2 /2 = c(0)/2, so 
that the parameter #1 is exactly a half-width of the main 
correlation peak. 



2.2 Basic relations 

It is known (see, e.g., iRoddierl Il98lft . that under the ap- 
proximation of weak perturbations, the scintillation index 
measured at the Earth's surface, is defined by the distri- 
bution of the structural coefficient C 2 (z) along the line of 
sight throughout the whole atmosphere using the following 
equation: 



= / Cl(z)Q(z)Az, 

J A 



(3) 



where 2 is the distance of the turbulent layer. For a measure- 
ment in the zenith, the distance equals to the layer altitude 
h above the site. Similarly one can define the cross-index 
(covariance): 



c(9) 



Cl(z)R(z,6)Az, 



(4) 



The Q(z) and R(z,6) are the weighting functions relating 
the OT intensity inside the layer to its effect in the observ- 
ing plane. Small-amplitudes approximation works always for 
large apertures, so this assumption is not restrictive for what 
follows. 

The technique f or computing the Q(z ) was described 
previously (see, e.g., iTokovinrnl l2002i . [2003) in details. For 
Kolmogorov model of isotropic and locally homogeneous tur- 
bulence with 2D spatial power spectrum of wave phase per- 
turbations oc f 11 ^ 3 (where / is the modulus of the spatial 
frequency), there is a known expression for measurement 
with zero exposure time: 



Q(z) = 9.6l/ r 11/3 S(z,f)A(f)fdf. 



(•») 



where A(f) is the aperture filter which is an Airy function 
for circular receiving apertures. The Fresnel filter S(z, f) 
describes the propagation of the amplitude perturbations 
of the wave front. In the case of monochromatic light with 
wavelength A, S(z, f) = sin 2 (n\zf 2 )/\ 2 . The scale factor of 
the filter is the Fresnel radius r-p = (Az) 1 / 2 . The integrand 
in this formula represents the spatial power spectrum of the 
scintillation. 

Following the Wiener-Khintchine theorem, the covari- 
ance of the scintillation between two points separated by 
a distance A is the result of the Fourier transform of the 
spatial spectrum. All the previous expressions were aver- 
aged over the polar angle <f> since the OT spectrum and 
all the spectral filters are axi-symmetric. Accordingly, to 
calculate the R(z,8), it is convenient to use the repre- 
sentation of the Fourier transform in polar coordinates: 
X(A) = 2tyJ X(f) f J (2ttA/) df, which immediately leads 
to the expression for the covariance weighting function: 

R(z, 6) = 9.61 r r 11/3 S(z, f) A(f) J (2n6zf) f df. (6) 
Jo 

where the distance A was replaced with z9 from simple ge- 
ometric considerations. 
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It can be noted that R(z, 0) = Q(z) and the normalised 
function of the angular correlation (ACF) in the layer is 
defined as p(z, 8) — R(z, 6)/Q(z). For the whole atmosphere, 
the observed ACF g(9) represents a normalised cross-index: 



q(6) 



<8± 

„2 



)Q{z)p{z,6) dz. 



(7) 



2.3 Large aperture specifics 



For large telescope with aperture of diameter fl > rp, the 
spectral filter S(f) entering in the integral ([6]) can be simpli- 
fied: the filter A(f) overlaps with S(f) a little because the 
first passes the frequencies below ~ 1/D and the second has 
main maximum near 1/rp. In the low frequency range, we 
may repla ce sin 2 (7rAz,f 2 ) by (rrAz/ 2 ) 2 without changing the 
integrand (Roddicr 1981). This technique is even more suit- 
able when calculating the covariance because the function 
Jo suppresses additionally high frequencies. 

This simplification leads to important results: 1) the 
integral dependence on the propagation distance z becomes 
oc z 2 , and 2) the dependence on light wavelength disappears 
altogether. The latter means that the scintillation becomes 
almost achromatic, therefore there is no need to use the 
exact expression of the Fresnel filter (|Tokovininll2003l l for 
polychromatic light. 

For subsequent calculations, we replace the modulus / 
by the dimensionless frequency q = fD. Then the approxi- 
mate Fresnel filter assumes the form of S(z, q) = n 2 z 2 D~ 4 q 4 
and the aperture filter changes to A(q) — [2 Ji(nq) /nq] 2 for 
clear circular aperture. After these transformations, the in- 
tegral (f5]) is easily evaluated and is equal to 0.4508, which 
leads to the usual formula for Q(z) = 17.34 D~ 7 ^z 2 . 

Most astronomical telescopes have a central obscura- 
tion caused by a secondary mirror, for an annular pupil, the 
aperture filter A(q) is as follows: 



A(q,e) 



Ji(nq) — eJijneq 



1 



(8) 



where e is equal to the ratio of inner to outer diameters. The 
use of this aperture filter affects only the numerical factor 
in the expression for Q(z). 

2.4 Short exposure regime 

Here we consider the angular correlation of the scintillation 
in the case bf of exposure time r so short that the OT carried 
by the wind w in the atmosphere does not have time to move 
across a noticeable distance. For small apertures, the short 
exposure condition is r <C rp/w (the order of few 0.001 s), 
while for large telescopes, r <C D/w (~ 0.01 or even ~ 0.1 s). 
The latter case is well described by the covariance function 
([6]) which becomes, after transition to frequency q: 



R(z,6) = 38.4 D 



-7/3 2 
' Z 



-2/3 



[J 1 (7rg)] 2 Jo(27r 7 q)dg, (9) 



where the integral depends only on the normalised angle 
7 = 9z/D. 

For a single layer at distance z, the ACF p(z, 6) depends 
only on 7 as well: 



p( 7 ) = 2.218 



-2/3 



[Ji(nq)] J (2nyq)dq, 



(10) 
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Figure 1. Angular correlation of the scintillation depending on 
the normalised angle 7 = dz/D for short exposure regime. The 
curves from top to bottom are computed for telescopes diameter 
of 1 m, 2.5 m and 6 m for turbulence altitude of 15 km (thin lines) 
and 5 km (thin dashed lines) are shown. The thick dashed line 
depicts the approximation mil l. In the box, a zoom of the ACFs 
at low angle is shown. 



and can be calculated analytically through generalised hy- 
pergeometric fun ctions (see equation f AID . A similar result 
was obtained in (Irlickson fc Lanzettal 120041 ) where the spa- 
tial correlation of the scintillation from extended object was 
considered. There the section of the cone of rays from the 
Moon played a role of the telescope aperture. In that paper 
non-normalised functions were calculated, which differ from 
our expressions by factor 0.4508. 

The result of the integration can be approximated by 
its series expansion: 



P(7) 



1.0967 



2/3 



0.243 7 + 0.231 7 



8/3 



(11) 



This approximation is shown in Fig. [T] on which exact func- 
tions, calculated by numerical integration of equation @ 
with A = 500 nm, are also plotted. One can see that, at 
8z < rF, the exact curves are above the approximation due 
to the simplified Fresnel filter in formula @. 

The value of the ACF at small angle is defined by the 
diffraction. It is easy to understand this by considering the 
original expression ((6} as the convolution of the ACF for 
an infinitely small aperture with the auto-correlation func- 
tion of the entrance pupil. The knowledge of the ACF at 
low angle is essential in problems where a strong correlation 
is expected. According to Fig. [TJ the diffraction effect for a 
2.5 m diameter telescope can increase the correlation coef- 
ficient, e.g., from 0.9 to 0.95 bf compared to the simplified 
expression, which corresponds to a two-fold reduction of the 
variance of the differential signal. 

However, the expansion approximates well the ACF in 
the domain 7 < 1. The angle at 50 per cent correlation is 
obtained from JTT} for one turbulent layer: 



h = 0.30 — . 

z 



(12) 



It is possible to estimate 9i by choosing the altitude of the 
turbulence which produces the main impact on the correla- 
tion. 
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Figure 2. Short exposure ACFs dependening on the normalised 
angle 7 = 9z/D for large telescopes with different central obscu- 
ration e. The thick line depicts e = 0, then e is equal to 0.1, 0.2, 
0.3, 0.4 and 0.5, respectively. 



Table 1. Relative dependence of the angular correlation radius 
01 on telescope central obscuration e. 



e 


0.1 


0.2 


0.25 


0.3 


0.35 


0.4 


0.45 


0.5 


9(e) 


0.95 


0.81 


0.72 


0.67 


0.62 


0.57 


0.52 


0.47 



To calculate the ACF g(9) related to the whole atmo- 
sphere, the formula is to be applied. Assuming Q(z) oc 
z 2 , we obtain that g(9) represents the average of the p(z,9) 
over all layers weighted by C 2 (z)z 2 . The angle 9\, denned 
as q(6\) — 1/2, is, as expected, directly proportional to the 
telescope diameter. 

It is well known that teles cope central o bscuration af- 
fects the scintillation power (|Yound 1 19671 ; iDravins et alj 
1998). It turns out that it acts upon the angular correlation 
as well. Fig. [2] shows the ACFs calculated numerically for 
different e. For large e, the shape of the curves becomes very 
complicated. These curves resemble the modulation transfer 
functions for annular apertures, due to the similarity of their 
respective mathematical expressions. 

The correlation angle for various central obscurations is 
obtained by multiplying 9i computed for circular aperture 
by some factor g(e) listed in Table [T] Note that when the 
inner diameter of the entrance pupil eD becomes comparable 
to the tf (telescope diameter < 0.5 m), the step on the 
curves p(-y) almost disappears. 



2.5 Long exposure regime 

In the previous section the averaging induced by the OT 
motion during the exposure was negleted. However, the ex- 
posure t is usually such that the wind can shift the OT over 
distances exceeding the telescope aperture size wt > D. 

The wind shear of the OT has been introduced in the 
context of the transition from the spatial correlation to 
the spatio-temporal correlation (jRocca. Roddier fc Verninl 
ll974l ; ICaccia. Azouit fc Verninlll987r i. With respect to spe- 
cific photometric issues, the dependence of the scintilla- 



tion on the exposure time was st udied in |Younsl 1 19671 : 
IDravins et al.lll998l ; lKornilovll2011b|). 

In a previous paper |Kornilovll2011lJ) we replaced tem- 
poral averaging over the exposure time with spatial averag- 
ing over an extended area calculated a s the convolution of 
the aperture D and the wind shear wt |Tokovinidl2"002r i. It 
is worth noting that in the case of a large wind shear, the 
"frozen turbulence" hypothesis may be not be fully valid 
due to the temporal evolution of the turbulence. This does 
not matter in the present case because the spatial spectrum 
remains unchanged. 

The angular correlation of the scintillation in a given 
direction depends on the relative orientation of the wind 
w(z). Indeed along the wind direction, the averaged region 
of the wave front has an area of about wtD. In the case of 
transverse displacement at distance A, the non-overlapping 
area « 2A/D is much larger than in the case of displacement 
along the wind w 2A/wt. 

Since there is a preferred direction, the problem ceases 
to be axi-symmetric and a two-dimensional integration over 
the spatial frequency is required. We carry out this inte- 
gration in polar coordinates introducing, as dimensionless 
variables, the frequency q = Df, aperture shift 7 = 9z/D, 
and wind shear ui = wr/D. We rewrite the expression 
in the form of a two-dimensional integral and we add the 
factor sinc 2 (cjg cos<^>), which describes the averaging by the 
wind along the a: axis l|Tokovininll2002l ; lKornilovll2011bl lah: 

/>00 

R(z, 9) = 9.6LD 5/3 / q~ 11/3 S{z, q)A(q) q dq 
Jo 

1 f 2w 2 

x — / cos(27r7q cos(<j!> — i/0)sinc (jnoq cos (j>) &(j>. (13) 
2ti" Jo 

Here tp is the position angle of the second object in this 
coordinate system. Note that the function sine x — sin x jx 
is not normalised. Let us denote £(g, w,7, tp) the integral 
over the polar angle <j>. 

The numerical evaluation of the expression (|13|) con- 
firms the strong dependency of the angular correlation to 
wind direction. However, as the wind can move in any di- 
rection relative to the pair of compared stars, the mean effect 
of wind smoothing is especially important. To evaluate it we 
may average the expression (|13[) over the angle ip without 
loss of significance. 

In the integrand of C(q,ui,^/,tj}), the multiplicand sine 
does not depend on the angle tp, so we can change the order 
of integration and take it out of the inner integral. The in- 
ner integral does not depend on <j> and equals 2nJo(2n'yq). 
Integration over (j> gives the spectral filter of win d smooth- 
ing 71 (u>q). An analytic representation is given in |Kornilovl 
l2011bl ). Finally we have 



£(q,u,7) =Ti(uq)J (2ir-yq). 



(14) 



Since, in the long exposure regime, the T\(u}q) has l/nuiq for 
asymptote, the average covariance function can be written 
in the form 



R(z, 



9.61 



£8/3 



-11/3 



S(z,q)A{q)J Q {2-K lq )Aq. (15) 



The equation can be simplified further using large aperture 
spe cifics (Sect. 12.31) and, for 7 = 0, turns into expression (7) 
in (|Kornilovl l2011b): Q(z) = 10.66 D~ A ' Z z 2 [wt)~ 1 . After 
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Figure 3. Long exposure (thick lines) versus short exposure func- 
tions (thin lines) ACFs averaged over all wind directions as a 
function of the normalised angle 7 = 6z/D. The clear circular 
aperture case is shown with solid lines, the annular aperture case 
e = 0.43 with dashed ones. 

normalisation, we derive the average ACF p(z,~/): 

r 00 

p{z,9) = 1.15 / q- 5/S [Ji(Trg)] 2 J (27r 7 g)dg. (16) 
Jo 

The integrand in (|16[1 differs from Sect. 12.41 in multiplicand 
q,-5/3 instead of q~ 2 ^ 3 . That reduces further the share of 
high-frequencies in the scintillation spectrum. 

We see that the function p(z, 7) depends only on the 
normalised angle 7 and does not depend on wind shear. Ac- 
cordingly, the correlation angle 9i is also independent on 
wind speed and exposure time as long as wt > D. 

The ACFs p(y) are shown in Fig. [3] in comparison with 
short exposures functions. The functions become smoother 
and the effect of the central obscuration becomes less notice- 
able. It is seen that the angle 8\ increases by a factor of 1.93 
(7 = 0.57) for clear aperture, and by a factor of 3.35 for a 
telescope with e = 0.43. It is particularly important that the 
behaviour near has radically changed. For the correlation 
of 0.9, the normalised angle is increased from 7 = 0.035 to 
7 = 0.178, which means an increase of the angular distance 
from 1 to 5 arcsec for a 2 m class telescope . These numbers 
are c onsistent with estimates obtained by l|Rvan fc Sandler! 
Il998h . 

The result of the integration of (|16|l is presented in Ap- 
pendixfX] equation (1A2[) . A good approximation of the ACF 
in the range 7 < 1 can be obtained using its series expansion: 

p{-y) w 1 - 2.79 7 5/3 + 1.78 -y 2 + 0.311 7 11/3 - 0.141 7 4 . (17) 

The previous equations are valid for clear circular aper- 
tures, but Fig. |3 shows that the introduction of a central 
obstruction does not change the value of the angle giving 
1/2 in the long exposure regime. 



3 PREDICTION OF THE ANGULAR 

CORRELATION OF THE SCINTILLATION 

We apply the above development to the vertical profiles 
of the OT monitored with the MASS/DIMM instrument 




10 20 30 40 

9, arcsec 

Figure 4. Example of the short exposure ACFs computed for real 
OT vertical profiles for the SAI 2.5 m telescope. 100 randomly 
chosen records are drawn, the bright dots denote the median of 
1684 individual records. 

ijKornilov et al.ll2003l . 120071 ) at Mt Sh atdzhatmaz (2127 m 
above sea level) in the North Caucasus (jKornilov et al.l2010l ) 
to estimate the angular correlation for the whole turbulent 
atmosphere. 

The database for the period November 2007 - October 
2009 consists of w 85 000 profiles. Each profile was obtained 
as a result of one-minute measurement and the restoration 
was performed on a 13-layers grid composed of the surface 
layer and 12 layers centered at 0.5, 0.71, 1.0, 1.41, 2.0, 2.82, 
4.0, 5.66, 8.0, 11.3, 16.0 and 22.6 km. The output of the 
MASS profile restoration is, for each altitude, the intensity 
AJi equal to the integral of the index of refraction structure 
coefficient within this layer: C^{zi)Sz. 

During restoration, the top node of the grid (22.6 km) 
accumulates the noise of the inversion process, and its out- 
put intensity may be somewhat overstated. For this reason, 
we decrease by half the intensity of this layer in further es- 
timates. The effect of such a correction does not exceed 10 
per cent. 

The results of the simulations should be used in plan- 
ning of future high-precision photometric observations and 
in the design of the SAI 2.5 m telescope equipment. That 
is exactly why in this paper we often use parameters of this 
particular telescope: entrance aperture diameter 2.5 m and 
central obscuration e = 0.43. 

3.1 Angular correlation in the short exposure 
regime for the SAI 2.5 m telescope 

In order to predict the angular correlation of the scintilla- 
tion on the SAI 2.5 m telescope with short exposures, the 
corresponding covariance function was calculated for each 
OT profile by the formula 

12 

c(e) = ^z*AJip(ezi/D), (18) 

i=0 

and then normalised by c(0) to obtain the ACF g(9) and 9i, 
the angle at which g(8i) = 1/2. We used the exact ACF in 
the layers p(9zi/D), calculated for each node Zi by the for- 
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Figure 5. Cumulative distributions of the correlation angle Q\, 
predicted for short exposures on the SAI 2.5 m telescope. The 
solid and dashed lines correspond respectively to the distributions 
for the telescope and for a clear circular 2.5 m aperture. 

mula (J6]), to take into account the central obscuration and 
diffraction effects. Fig. [4] shows randomly selected resulting 
ACFs. The median ACF for enlarged sample is plotted as 
well, showing the characteristic step due to central obscura- 
tion. Also note the long (30 — 40 arcsec) pedestal attributed 
to the presence of low-altitude turbulence. 

The median of the 6\ distribution on Fig. [5] is equal to 

7.1 arcsec. The first and third quartiles of the distribution 
are equal to 6.4 and 8.0 arcsec, respectively. Replacing 9i/D 
ratio in formula <|12[) . we estimate the typical effective alti- 
tude of the OT, after correction for the central obscuration, 
to w 12 km. 

3.2 Angular correlation in the long exposure 
regime for the SAI 2.5 m telescope 

Although, the normalised averaged ACF does not depend 
on wind, the weight of each turbulent layer does because 
the scintillation power is proportional to (wr)^ 1 for long 
exposures. Therefore, in addition to the OT profiles, some 
information is required on the vertical distribution of wind 
speed. For the prediction, we used the vertical wind pro- 
file \V (zj)} provided by t he NCEP/NCAR meteorological 
model (jKalnav et al.lll996r ) averaged during the observation 
period. 

The "observed" covariance function was constructed us- 
ing the formula 

12 . 

c (0)=X>?^rP(0V-D): (19) 

i— 

where p(6zi / D) in each layer was calculated using the exact 
formula (|15(l . Then, as in the previous case, the function was 
normalised and the angle 8i was determined. 

The cumulative distribution of 9\ for all OT profiles 
is shown in Fig. [6] Compared to the short exposures case, 
the correlation angle is significantly larger with a median 
of 20.3 arcsec, and the first and third quartiles of 19 and 
23 arcsec, respectively. The distribution is very narrow, and 
these numbers can be used with confidence for planning the 
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Figure 6. Cumulative distributions of the correlation angle 9\, 
predicted for long exposures on the SAI 2.5 m telescope. The 
dashed line shows the case without considering wind profile. 




Figure 7. Example of the long exposure ACFs computed by real 
OT vertical profiles for the SAI 2.5 m telescope. 100 randomly 
chosen records are shown, the bright dots denote the median of 
1618 individual records. 



observations. Neglecting the wind distribution changes the 
correlation angle by w 10 per cent only. 

Examples of individual ACFs and their median function 
are shown in Fig. [7] As in the previous section, the median 
ACF is computed as median values of individual functions 
for fixed 8 and respective angle Q\ is computed for q(9) = 
1/2. This procedure leads to a slightly different result from 
the simple median of all the individual 6\ . 

In general, the median ACF is similar to the function 
for a single layer (Fig. [3} ■ Features near 6 — are defined by 
the high and quite stable turbulence. A strong correlation 
of 0.9 is observed for objects at an angular distance of w 5 
arcsec. When 6 « 10 arcsec, the correlation coefficient can 
be about 0.75, which corresponds to twofold reduction of the 
scintillation noise with respect to the distance equal 6\ and 
fourfold compared to the case of far apart objects. 

For this photometric regime, the effective altitude h e g = 
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0.57 D/9\ is 14 km, which corresponds to a barometric alti- 
tude of 16 km. 



4 ANALYSIS OF THE SCINTILLATION NOISE 
REDUCTION METHODS 

4.1 Method of conjugated photometry 

IOsborn"eT al, (2010) proposed a method for differential pho- 
tometry where the entrance pupil of the telescope is conju- 
gated to the most intense turbulent layer. According to the 
authors, this operation increases the correlation of the scin- 
tillation and consequently reduces the scintillation noise in 
the differential signal. Such a geometry of light propagation 
increases the contribution of the low-altitude turbulence to 
the scintillation but the gain in the correlation is expected 
to be larger. 

However real OT vertical distributions may be substan- 
tially different from the simplified model used in the paper. 
First of all, the maximum contribution in the scintillation 
is not always produced by the most intense layer. Because 
of the z 2 weighting factor, higher layers with weaker OT 
can dominate. Secondly, in the case of short exposures, it 
was shown that the correlation is reduced very quickly even 
for small angular separation. This means that even a weak 
high-altitude turbulence can limit strongly correlated scin- 
tillation to a very narrow region of the sky. For the case of 
long exposures the ACF peak is broader but a role of slowly 
moving OT in the surface layer grows up. In any case, an 
amplification of the scintillation, caused by the transfer of 
the entrance pupil to the conjugate altitude z* , is very sig- 
nificant. 

To analyse the situation, we define 8 C as the angle at 
which bf conjugated photometry decreases the scintillation 
noise power by a factor of 2 compared to conventional dif- 
ferential photometry (assuming uncorrelated scintillation): 

s 2 = 2a* 2 (l- *(8 c )) or g*(0 c ) = l _ _fL, (20) 

where s 2 is the scintillation index under a given state of 
the atmosphere, a* 2 and Q*{8) are the scintillation power 
and the ACF for the whole atmosphere in the case of the 
conjugated pupil. Using the e xpression for the scintillation in 
a case of virtual propagation l|Fuchs. Tallon k, Verninlll998l ) 
and our data on the OT vertical distribution, we calculated 
the correlation as 

1 12 

Q*{6) = ]T AJi(zi - z*fp{0zi/D). (21) 

i=0 

Evidently, for objects closer than 8 C , the gain in noise sup- 
pression is greater than 2. Computing p(0Zi/D) we did not 
take into account the loss of a part of the telescope aperture 
in conjugated pupil mode. Increasing the central obscura- 
tion additionally amplifies the scintillation and reduces the 
correlation. 

Distributions of these estimates are presented in Fig [5] 
for short exposures. Fixed z* distributions show that in our 
conditions (when the surface layer contains w 0.65 of total 
turbulence), the conjugation to 5 km altitude is more prof- 
itable because it minimises the scintillation from the surface 
turbulence. Thus, when z* = 5 km, the median 8 C reaches 




Figure 8. Distribution of the angles 9 C calculated with real OT 
profiles for different conjugation altitudes z* (thin lines) in km. 
The thick solid line depicts the adaptive selection of the z* min- 
imising the scintillation in the free atmosphere, the dash-dotted 
line corresponds to the adaptive selection of the most intensive 
layer. The thick dashed line is the same as solid, but for the scin- 
tillation noise reduction by 3 times. The vertical line corresponds 
to 0i. 

190 arcsec, and in a quarter of the time 8 C > 20 arcmin. A 
similar picture is observed in the case of an adaptive choice 
of z* , minimising the scintillation in the free atmosphere 
(z ^ 1 km). Results are worse when conjugating to the most 
intense free atmosphere layer. 

The distributions show the "toggle" character of the 
process: in the case of strong surface turbulence, the angles 
8 C are small and the application of the method is mean- 
ingless. But above a certain level, the technique becomes 
effective. An attempt to reduce the scintillation noise more 
(by 3 times) is shown by the dashed line: an adaptive choice 
z* provides the median angle w 35 arcsec only. 

In the long exposures case the main characteristics of 
the distributions are preserved but the correlation angles 
become smaller. E.g., for the adaptive choice of the conju- 
gation altitude, the median 8 C is « 40 arcsec, and 8 C > 2 
arcmin occurs only in 18 per cent of the time. The main rea- 
son for the deterioration in conjugated mode is an increase 
of the scintillation due to slow low-altitude turbulent layers. 

For a successful application of the discussed technique, 
data on the current OT and wind above the telescope are es- 
sential since these characteristics are very unstable. A study 
of the temporal evolution of 8 C shows that it can change 
by several times over the time scale of tens of minutes. 
Note, that the altitude **, providing the minimum amplifi- 
cation of the scintillation, coincides with the effective alti- 
tude of the turbulence in free atmosphere. The softw are of 
the MASS/DIMM instrument l|Tokovinin et alj|2003t) deliv- 
ers this parameter in real time. 

4.2 Ensembles 

Differential CCD photometry using a close group of stars, 
located around the object of interest, was proposed primar- 
ily to minimise the impact of variations of the atmospheric 
transparency in order to obtain high-precision (better 0.001 
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magnitude) data for the study of stellar mic ro-variability 
l|Brown fc Gillilandl Il994l ; iHeaslev et all Il996l ). This tech- 
nique was a lso found useful for observation of systems with 
exoplanets (lEverett fc Howellll200ll ; ISouthworth et al.l l 2009l ; 
M ann, Gaidos & Alderindl201ll V 

The ensemble of the comparison stars should be fairly 
compact as the spatial correlation of the fluctuations in the 
extinction caused by, e.g., cirrus, has a radius no more than 
a few arcminutes. This issue is poorly understood, but ev- 
ery astronomer can easily verify this by looking at clouds 
patterns passing over the Moon. 

It was also noted that the ensemble reduces the contri- 
bution of the stellar scin tillation to the error budget of the 
measurements (see, e.g.. iGilliland et al.lll993h . We investi- 
gate in this section the potential of this technique invoking 
the angular correlation of the scintillation. First we discuss 
the procedure for comparing the measured signal. 

There are two obvious ways to compare the flux from 
the object 7, to the total flux Is of N stars in the ensemble: 

I - 
o=Y' where J s = 22 Ij, (22) 

i 

and/or alternatively, in the space of magnitudes 

1 N 

Am = m, - — 2~2 m i ( 23 ) 

3 

In general, to each star in the ensemble can be attributed a 
certain weight. Evidently, the formula (|22[) shows the con- 
tribution of multiplicative noise for fainter stars with less 
weight. 

We write the analogue of the formula ((1]) for the ensem- 
ble: 

2 

2| = s 2 [l - 22~2w j g(e*, j ) + 2~2 w2 j + 2 2~2 w i w kQ(hk)] , 

j j h,j<k 

(24) 

where Wj = Ij/ It, is the relative weight of j-ih star in the 
fluxes. In the magnitudes space, all the stars in the ensemble 
are equally accurate, and 

3 k,j<k 

In these expressions, the sum in brackets can be interpreted 
as an amplification of the scintillation noise. The second 
term in the sum describes the correlation of the star un- 
der study with the comparison stars, the third and fourth 
term are the variance of the mean brightness of correlated 
comparison stars. 

For clarity, we are going to analyse the latter formula. 
If the distances between all the stars are much larger than 
8i, the resulting variance is reduced to s 2 (l + JV -1 ), i.e. an 
increase of the number of the stars reduces the scintillation 
noise. However, if close pairs appear in the ensemble, further 
improvement is slowed down by the fourth term. Nonetheless 
in the limit, when the ensemble covers the sky area much 
larger than 0\, it can be regarded as an extended object, 
which has no scintillation. 

When the ensemble of stars are at distances from the 
object comparable to the radius of the angular correlation, 
the second term, the average correlation coefficient, is to be 
taken into account. Winning situation is achieved when the 



Table 2 . Amplification of the scintillation noise for different con- 
figurations of star ensemble. See text for explanation. 



R 


N = 1 


2 


3 


4 


6 


8 


oo 


0.5 6»i 


0.57 


0.32 


0.27 


0.26 


0.24 


0.24 


0.24 


1^1 


1.00 


0.63 


0.53 


0.50 


0.48 


0.47 


0.46 


2 6>i 


1.49 


1.05 


0.91 


0.85 


0.80 


0.79 


0.77 


4 0i 


1.78 


1.30 


1.15 


1.08 


1.01 


0.98 


0.93 



object under study is located in the centre of ensemble and 
is surrounded by a large number of comparison stars. 

As a numerical example, let us consider the ensemble 
of TV stars arranged in a circle of radius R around the star 
under study. The amplification for different configurations 
is given in Table [2] The median curve from Section 13.21 is 
used as the ACF with the angle 0\ = 23 arcsec. 

The first row of the table reflects the situation of a 
crowded stellar field when the typical distance between the 
stars (12 arcsec in our example) is less than the angle of the 
correlation. We see that in this case there is no reason to 
have more than 4 comparison stars. Further increase in their 
numbers do not improve the situation even if they are close 
to the star. Adding new stars at larger distance, we reduce 
the second term in (|25[) and increase the amplification. 

The last line describes the case of rarefied stellar ensem- 
ble, the typical distance between stars is 1.5 arcmin. When 
using a single comparison star, the amplification factor is 
« 1.8. However, we can reduce its value up to 1, correspond- 
ing to a differential measurement at distance 9\ , by choosing 
6 to 8 comparison stars. 

Calculations of the scintillation power from OT profiles, 
give the median scintillation noise equal to 250 /xmag for 
60 s integration on the SAI 2.5 m telescope. It means that 
the scintillation noise can be as small as 120 fim&g using 
the ensemble technique, to be compared to 350 ^imag with 
trivial differential photometry. 



5 DISCUSSION 

All the equations in the paper are written for the distance 
z to turbulent layer and not for its altitude h above the 
observatory. The dependency of the scintillation power on 
the zenith distance Z (or air mass AIz = secZ) is obtained 
by replacing z at h sec Z in the formulae ([3]) and Q what 
leads to oc M|. In the long exposure regime, the ex ponent 
varie s from 3 to 4, dep ending on the wind direction (|Yound 
1 19671 : iKornilovllinTrJ ). 

In any case, this is a common factor for all the layers 
and it does not change their relative contribution to the scin- 
tillation. Therefore, to identify the dependency of 0\ on air 
mass, we can use the formula (|12|l which implies 6\ oc M^ 1 . 
Thus, for differential photometry far from the zenith, not 
only the scintillation power greatly increases, but also the 
correlation of the scintillation decreases. In the two limit- 
ing cases, when the sources are either very close (q « 1) or 
very far (q w 0), the variance of the differential signal in- 
creases together with the s 2 . However, for distances of order 
0i (q « 0.5), o\ m grows faster. 

We did not touch the absolute intensity of the scintilla- 
tion, analysing the comparative effects, all other conditions 
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being equal. The atmosphere above the tropopause is gen- 
erally stable enough, however a growth of s 2 occurs in pres- 
ence of strong turbulent layers at the tropopause altitude 
(jet streams) or below in form of wind wake from mountain 
ridges or severe weather phenomena. In such cases, the angle 
of the correlation should also increase. Similar dependence 
is observed in our data, but it is very blurred and Qi changes 
only at ~ 50 per cent when the s 2 varies by a few times. 

Since the typical correlation angle corresponds to a few 
meters in the high atmosphere, a deviation of the spatial 
spectrum of the refractive index from the Kolmogorov model 
is to be taken into account. It is known that the finite outer 
scale of the non-Kolmogorov turbulence greatly reduces the 
spatial spectrum already in a m eter range, what affe cts a 
size of image at large telescopes (|Martinez et alj|2010t ). 

In our case, this factor suppresses low-frequencies scin- 
tillation what is to lead to a shortening of the correlation 
distance. This effect resembles the central obscuration effect 
which amplifies the high-frequencies. The outer scale at high 
altitu des is large, with a characteristic length about 20 to 
50 m (|Dali Ali et alj|2010h . This gives reason to assert that 
the obtained results are applicable for telescopes 2 - 4 m 
class with confidence. 

The comparison of photometric methods to reduce 
the scintillation noise: ensemble of comparison stars and 
conjugate-plane photometry, leads to the following conclu- 
sions. The first method requires a close group of stars which 
can not always be found. But if there is a compact group of 
4-5 stars, a reduction of the scintillation noise is possible 
up to 3 - 4 times. The method is always applicable because 
is associated to the high-altitude turbulence. With larger 
telescopes, more extended sky area can be used. 

The feasibility of conjugate photometry depends on the 
turbulence in the lower atmosphere. Usually it is distin- 
guished by high instability. Our simulations show that in 
the sites where low-altitude OT normally dominates, three- 
fold gain in signal-to-noise ratio is not always attainable. 
But in certain atmospheric conditions, the comparison star 
can be chosen far from the target star (2 — 4 arcmin). The 
question is how to predict and detect these conditions? 



diet the angular correlation of the scintillation for photomet- 
ric observations at this observatory. All hereinafter referred 
numerical results relate to this telescope and this site. 

Since in case of a large telescope, the ACF depends 
only on one parameter: the normalised angle 7 = 6z/D, 
it is possible to scale the obtained estimates for telescope 
with different diameter. Central obscuration in a telescope 
affects the correlation of the scintillation for short exposures 
but not for long exposures. 

Simulation for short exposure regime (fast photometry) 
predicts the median of the 9i as large as 7 arcsec. In the long 
exposure regime, the ACF averaged over all wind directions, 
does not depend on wind speed and exposure. In this case, 
estimated median of the correlation angles 81 amounts to 
m 20 arcsec. 

Study of the conjugate-plane photometry shows that 
the method depends critically on the intensity of low- 
altitude OT. When the lower atmosphere is calm, one can 
use the comparison star at distance greater than several ar- 
cminutes. Analysis of the ensemble photometry shows that 
the compact, but quite resolved, ensemble of four compari- 
son stars is enough to get the reduction in the scintillation 
noise power by ~ 8 times. Photometric observation of the 
different ensemble configurations could provide the neces- 
sary experimental verification of the proposed formulae and 
the performed simulations. 
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6 CONCLUSION 

In this paper, we discussed the problem of correlation of 
the stellar scintillation from two point-like sources, closely 
spaced in the sky. Using the standard model of the optical 
turbulence (Kolmogorov spatial spectrum of refractive in- 
dex fluctuations) in the approximation of the weak pertur- 
bations, two situations were considered for large telescope 
aperture (much larger than the Fresnel raduis tf): 1) ob- 
servations with short exposures, 2) observations with long 
exposures (typical photometric practice). 

Analytic expressions obtained for the ACFs, have 
lengthy description, therefore a simple power-law approxi- 
mation, well describing the functions, has been proposed. 
As the measure of the angular correlation, the angle 61, at 
which the correlation coefficient is equal 1/2, was adopted. 

Based on previously obtained data on the vertical distri- 
bution of the OT in the atmosphere above Mt Shatdzhatmaz 
in the Northern Caucasus, the simulation experiments have 
been performed for the SAI 2.5 m telescope in order to pre- 
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APPENDIX A: EXACT FORMULAE FOR ACFS 
Al Short exposures 

In the case fl > tf (irrelevance of diffraction effect) with- 
out central obscuration in the short exposure regime, the 
integration of the (j 10|) gives: 

/9 (7) = 3F 2 (-f,|,|;| ) l; 7 2 )- 

2/3 rp I 1 1 3.4 4 . 



£i7 2/ %M-U,!;i,i; 7 2 ), (T< 1 ) 

'7 7 3 . . 

.6' 6' 2' - 



p(7)=^27- 7/3 3iMU,§;2,3; 7 - 2 ), (7>1) (Al) 



where 



Di = 15y|2^ r (|) p4 (|) = L0965j 

D 2 = -|ff = -0.0259. 

288 v rr 



A2 Long exposures 

In the case of large telescope without central obscuration in 
the long exposure regime, the integration of the (|16|) results: 

„( 7 ) = 3 F 2 (-1,-1, |; i 1; 7 2 )- 

-£ l 7 5/3 3f 2 (-ii|;f,f;7 2 ), (7<i) 

p(7) = £2 7~ 4/3 3f2(|,|,|;2,3; 7 - 2 ), (7>1) (A2) 



